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ABSTRACT 


The  expansions  of  algebraic  functions  can  be  computed  "fast"  using  the 
Newton  Polygon  Process  and  any  "normal"  iteration.  Let  M(j)  be  the  number 
of  operations  sufficient  to  multiply  two  jth  degree  polynomials.  It  is 
shown  that  the  first  N terms  of  an  expansion  of  any  algebraic  function  de- 
fined by  an  nth  degree  polynomial  can  be  computed  in  0(n(M(N))  operations, 
while  the  classical  method  needs  0(Nn)  operations.  Among  the  numerous  ap- 
plications of  algebraic  functions  are  symbolic  mathematics  and  combinatorial 
analysis.  Reversion,  reciprocation,  and  nth  root  of  a polynomial  are  all 
special  cases  of  algebraic  functions. 
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1.  INTRODUCTION 


Let 

(1.1)  P(W,z)  = An(z)Wn  + ...  + A0(z), 

where  the  A^(z)  are  polynomials  over  a field  A.  In  general  we  shall  take  A 
to  be  the  field  of  complex  numbers;  an  exception  being  Section  7.  (Many  of 
the  results  hold  for  an  algebraically  closed  field  of  characteristic  0.) 
Without  loss  of  generality  we  assume  Aq(z)  » 0 and  A^(z)  j=-  0.  Capital  letters 
will  denote  polynomials  or  series;  lower  case  letters  will  denote  scalars. 

The  zero  of  (1.1),  a function  S(z)  such  that  P(S(z),z)  = 0,  is  called 
the  algebraic  function  corresponding  to  P(W,z).  Let  z ^ be  an  arbitrary 
complex  number,  finite  or  infinite.  It  is  known  from  the  general  theory  of 
algebraic  functions  that  S(z)  has  n fractional  power  series  expansions  around 
z^.  by  the  computation  of  an  algebraic  function  we  shall  mean  the  computa- 
tion of  the  first  N coefficients  (including  zero  coefficients)  of  one  of  its 
expansions.  (This  will  be  made  precise  in  Section  3.)  The  problem  we  study 
in  this  paper  is  the  computation  of  one  expansion  of  the  algebraic  function. 
Our  results  can  easily  be  modified  for  computing  more  than  one  expansion  or 
all  expansions  of  the  algebraic  function. 

As  described  in  most  texts,  the  classical  method  computes  algebraic 
functions  by  comparison  ol  coefficients.  It  is  not  difficult  to  show  that 
the  method  can  take  0(NU)  operations,  where  n is  the  degree  of  P(W,z)  with 

respect  to  W.  Hence  the  classical  method  is  very  slow  when  n is  large. 

The  main  result  of  this  paper  is  that  every  algebraic  function  can  be 
computed  fast.  Let  M(N')  denote  the  number  of  operations  sufficient  to  multiply 
two  Nth  degree  polynomials  over  the  field  A.  Let  C(N)  be  the  number  of  opera- 
tions needed  to  compute  any  algebraic  function.  We  prove  that 


C (N)  = 0(nM(N) ) . 


Since  M(N)  = 0(N2)  (or  M(N)  = 0(N  log  N)  if  the  FFT  is  used) , our  algorithms 
are  considerably  faster  than  the  classical  method  even  for  moderate  n.  It  is 
an  open  problem  whether  or  not  a general  algebraic  function  can  be  computed 
in  less  chan  0(M(N))  operations- 

The  "fast  computation"  of  the  title  is  because  the  coefficients  of  a 
"regular"  problem  can  always  be  computed  fast  by  iteration  (Section  5)  and 
the  general  problem  can  b.;  reduced  to  a regular  problem  (Section  6)  with 
cost  independent  of  N. 

Brent  and  Rung  [1976]  showed  that  the  cost  for  reversion  of  a polynomial, 

which  is  a very  special  case  of  ^n  algebraic  function  (see  discussion  later 

2 

in  this  section),  is  0((N  log  N)  M(N)).  We  stated  above  that  the  cost  of 
expanding  an  algebraic  function  is  0(nM(N)).  These  results  are  reconciled 
by  the  observation  that  we  are  considering  the  case  that  the  degree  n of 
P(W,z)  with  respect  to  W is  fixed  and  independent  of  N,  while  Brent  and  Rung 
considered  the  case  where  n = N. 

There  are  known  examples  of  fast  computation  using  Newton-like  iteration 
in  settings  such  as  algebraic  number  theory  (Bachman  [1964]),  power  series 
computation  (Rung  [1974],  Brent  and  Rung  [1976]),  and  the  Zassenhaus  construc- 
tion in  p-adic  analysis  (Yun  [1976]).  Fast  computation  of  algebraic  functions 
raises  certain  issues  not  present  in  these  other  settings;  see  especially 
Section  6.  As  we  will  see  in  Section  5,  there  is  nothing  special  about  Newton-like 
iteration;  any  "normal  iteration"  can  be  used. 

Although  the  complexity  results  are  stated  asymptotically,  Theorems  5.1 
and  6.1  give  non- asymptotic  analyses  of  the  algorithms.  Hence  various  non- 
asymptotic  analyses  can  also  be  carried  out. 

We  are  interested  in  the  computation  of  algebraic  functions  for  a num- 


ber of  reasons.  These  include 


1. 


A number  of  problems  where  fast  algorithms  are  known  are  special 
cases  of  algebraic  functions.  (More  details  are  given  below.) 

2.  There  are  numerous  applications.  For  example,  many  generating 
functions  of  combinatorial  analysis  and  functions  arising  in 
mathematical  physics  are  algebraic  functions.  The  integrands  of 
elliptic  and  more  generally  Abelian  integrals  are  algebraic  func- 
tions. See  Section  9 for  an  example. 

3.  Algorithms  for  expanding  algebraic  functions  are  needed  in  systems 
for  symbolic  mathematics  such  as  MACSYMA  (Moses  [1974]). 

4.  Algebraic  functions  are  of  theoretical  interest  in  many  areas  of 
mathematics.  These  include  integration  in  finite  terms  (Ritt  [1948]), 
theory  of  plane  curves  (Walker  [1950]),  elliptic  function  theory 
(Briot  and  Bouquet  [1859]),  complex  analysis  (Ahlfors  [1966],  Saks 
and  Zygmund  [1971]),  and  algebraic  geometry  (Lefschetz  [1953]). 
Algebraic  function  theory  is  a major  subject  in  its  own  right.  See, 
for  example,  Bliss  [1933]  and  Eichler  [1966], 

We  exhibit  special  cases  of  algebraic  functions  where  fast  algorithms 
are  known. 

A . Reciprocal  of  a polynomial : 

P(W,z)  « AL(z)W  - 1.  (See  Rung  [1974].) 

(Actually  Rung  uses  P(W,z)  = w”1  - A^z)  which  is  not  of  the  form 
(1.1)>  and  allows  A^(z)  to  be  a power  series.) 

B.  nth  root  of  a polynomial: 

P(W,z)  * Wn  - AQ(z) . (See  Brent  [1976  , Section  13]  where  the  ^(z) 
is  allowed  to  be  a power  series.) 
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C.  Reversion  of  a polynomial: 

Let  f be  a given  polynomial  with  zero  constant  term.  We  seek  a function  g 
such  that  f(g(z))  ■ z.  To  see  this  is  a special  case  of  an  algebraic  func- • 
tion,  let  f(x)  = anxR  + an  ^xn  * + ...  + a^x.  Then  we  seek  g(x) 
such  that  angn(z)  + ...  + a^g(z)  - z = 0.  This  is  an  instance  of 
our  general  problem  with  A^(z)  ® a^,,  i=l,...,n,  Aq(z)  = -z. 

See  Brent  and  Rung  [1976]. 


We  summarize  the  results  of  this  paper.  In  Section  2 we  show  that 
without  loss  of  generality  we  can  take  z^  = 0 and  assume  An(0)  f 0.  Nota- 
tion is  established  and  a few  basic  facts  from  algebraic  function  theory 
are  summarized  in  Section  3.  The  concept  of  normal  iteration  is  introduced 
in  Section  4 and  convergence  of  normal  iterations  for  regular  problems  is 
established  in  Section  5.  In  Section  6 we  state  and  analyze  the  Newton 
Polygon  Process,  which  reduces  the  general  problem  to  a regular  problem. 

A symbolic  mode  of  computation  with  exact  arithmetic 

is  introduced  in  Section  7.  Section  8 shows  that  C(N)  = 0(nM(N)).  In  Section 
9 we  give  a number  of  examples,  several  of  which  are  more  general  than  the 
theory  of  the  preceding  sections.  The  final  section  discusses 
extensions  of  the  work  presented  here. 

In  this  paper  we  analyze  algorithms  under  the  assumption  that  the  coef- 
ficient of  power  series  are  "non -growing",  e.g.,  all  coefficient  computations 
are  done  in  a finite  field  or  in  finite-precision  floating-point  arithmetic. 

An  analysis  dealing  with  variable- precis  ion  coefficients  is  yet  to  be  performed. 
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2.  PKELMIHAllY  TRAliSFORMATIOHS 


,-:o 


Recall  that  we  wish  to  compute  doe  of  the  expansions  around  zn  of  the 


V 


algebraic  function.  S (a)  corresponding  to 


V*  i 

A 


P(V,z)  - An(r)W“  + ...  t AQ(z>, 


;V~ 

■-  \;> 


‘A  sZ> 


f y4 


i.ei,  P(S(z),z)  = 0. ; In  this  section  we  show  that  *ffc«  tmo  sispU  tttts- 

■ ■ ■ i. » , * ; • ■ . , . ■* 

formations  we  need  only  deal  with  idle  case  that  s^  ■ 0 and  An(0)  0.  If  we 

/ - I 

transform  P(W,z)  to  P(W,z),  then  S(z)  is  defined  by  P(S(s),z)  * 0. 


y 


Consider  first  the  case  z.  » *.  Let 


(2.1)  P(W,z)  - zmP(W,|) 


where  m » max  (deg  A.) . By  definition,  an  expansion  of  S(z)  around  z.  * 
0 si  sh  ■ 


is  an  expansion  of 


around  zQ  ■ 0. 


Consider  next  the  case  that  Zq  is  any  finite  complex  number  * Define 


P(W:,z)  - P(W,z+zQ). 


o 


Oi 


6 


An  expansion  of  S(z)  aroum- ; the  origin  is  Ah  expansion  of  S(z)  around  i * *q. 
For  the  remainder  of  this  paper  we  shall  therefore  take  z^  ■ 0,  v- 


Let  An(0)  • 0.  Then  the  algebrcic  function  S(z)  corresponding  tc  P(W,z) 


has  one  or  more  expansions  with  negative  powers.  Using  the  following  trans- 
formation, we  need  only  deal  with  expansions  with  non- negative  powers.  It 
is  convenient  to  use  ord  notation. 


Definition  2 .1.  Let  A(z)  be  an  integral  or  fractional  power  series,  if  A(t)  ^ 0, 
then  ord (A)  denotes  the  degree  of  the  lowest  degree  term  in  A(z) . If  A(z)  - 6, 
then  ord (A)  ■ ®. 


Choose  non-negative  integers  u and  X to  satisfy  the  following  conditions: 


U + ord(An)  = nX, 

u+  ord(A^)  £ iX,  i=l,...,n-l. 


P(W  ,z)  = z^P(W/z\z). 

Then  the  coefficients  of  P(W,z),  A.(z),  are  polynomials  with  A (0)  0,  and 

i n 

5(z)  has  only  expansions  with  non-negative  powers.  Since  the  expansions  of 
S(z)  are  those  of  S(z)  divided  by  z\  it  suffices  to  compute  expansions  of 
S(z).  For  the  remainder  of  this  paper,  we  therefore  assume  tnat  A^CO)  * 0. 
(One  should  note,  however,  that  the  results  of  Section  5 hold  without  the 
assumption.) 


N 


A 


f 


wiSPWiCTJ  w mfm&T&zz  gg^  wBr^gw^g^^ 
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3.  FACTS  FROM  ALGEBRAIC  FUNCTION  THEORY 


We  introduce  some  notation  and  state  a basic  result  of  algebraic  func- 
tion theory  which  characterizes  the  expansions  of  the  algebraic  function 
corresponding  to 

P(W,z)  = An(z)Wn  + ...  + Aq(z). 


There  exist  r positive  integers  d^,...,dr  such  that  d1  + 
expansions  of  the  algebraic  function  are  given  by 


(3.1) 


e 

L,  ' i , 0 *i 

0=0 


+ d - n and  the 
r 


for  i=l,...,r  and  j=0, . . . ,d 1 , where  ^is  a primitive  d^th  root  of  unity 

and  the  s,  . are  complex  numbers.  For  each  L,  the  expansions  S , , 

i,  Hi  J-»j 

j=C, . . . ,d^-l,  are  said  to  constitute  a cycle. 

The  problem  considered  in  this  piper  is  to  compute  one  expansion  of  an 
algebra!1'  function.  For  notational  convenience,  let  the  expansion  be  denoted 
by 

CO  l 

S (z)  = ) S „ • zd. 

t I AJ 

Hence  our  problem  can  be  formulated  as  that  of  computing  the  value  of  d and 
the  coefficients  sQ,s  ,...  . (In  this  paper  S(z)  represents  either  an  alge- 
braic function  or  one  of  its  expansions,  depending  upon  the 
context.)  Note  that  since 


P(S(z) ,z)  = 0, 


we  have 


P(s0  > 0)  = 0. 


•8- 


Thus,  Sq  is  a zero  of  the  numerical  polynomial  P(W,0).  We  say  our 

problem  is  regular  with  respect  to  sQ  if  sQ  is  a simple  zero  of  P(W,0). 

(In  this  definition,  we  allow  A (0)  to  be  0 .)  For  a regular  problem,  we 

n 

have  d=  1,  that  is,  the  expansion  S(z)  is  an  integral  power  series.  In 
Section  5,  we  shall  show  that  a regular  problem  can  always  be  solved  by 
iteration.  In  Section  6,  we  shall  show  how  the  general  problem  can  be 
transformed  to  a regular  problem. 


-a.-'-.  ■■•*«.-  !.n  (p,  i.  <! ■•ul"«„'. w^, 


-)- 


A . NORMAL  ITERATIONS 

We  introduce  the  concept  of  a normal  numerical  iteration.  We  give  a 
novel  definition  of  the  order  of  a normal  iteration  which  is  convenient  for 
the  application  to  power  series  iteration.  In  the  following  section  we  will 
show  that  a normal  iteration  with  order  greater  than  unity  will  always  con- 
verge il  used  for  a regular  problem. 

Let  p(w)  be  the  numerical  polynomial  P(W,D),  let  s be  a zero  of  p(w),  and 
/ * \ * • \ 

let  e'  = w'  - s denote  the  error  of  the  ith  iterate.  To  motivate  the 
definition  of  normal  iteration  we  first  consider  two  examples. 

Example  A . 1 . Newton  Iteration 

W(i+D  = wd)  . , 

P,(w(i)'> 

From  the  Taylor  series  expansions 


/ (i>  \ . , , \ ( t)  . i>"(s  > , ( L > 2 
P(w  ) - p (s)  e + ~ — (e  ! + 


and 


|)’  (W<L)  - p’  (s)  + p"(s) 


U > 


we  have 

<4-l>  ^<i+1)  " vfe(<'u>)2+  ...  cj  ’ 

j=3 

where  the  c^  are  rational  expressions  of  the  derivatives  of  p at  s,  with 
powers  of  p1 (s)  as  the  denominators. 

Example  A. 2.  Secant  Iteration 


(i)  a-i) 
W(i+D  = w(0  _ „-w 


, (i) . , (i-  i 1 . 

p(w  )~i>(  w > 


> > i 

• P(w 


Using  the  Taylor  series  expansions  of  p(w‘  ')  and  p(wv  '),  we  obtain 


(4.2)  e(tfl)  - e<1)e(t"1)  -f 

K ’ 2p'(s)  6 + 


L cjJ t 
j+£*3 


where  the  are  rational  expressions  of  the  derivatives  of  p at  s,  with 
powers  of  p' (s)  as  the  denominators. 

Consider  now  a general  iteration 


..  (i+1)  , , (i)  (i-1)  (i-m). 

(4.3)  w'  = <!i(w  ,wv  , ...,w^  ) , 


which  is  defined  in  terms  of  rational  expressions  of  p and  its  derivatives. 
Assume  that  by  using  Taylor  series  expansions,  we  can  derive 


(4.4)  e(i+1)  - c. 

^ J0’***,Jm 

j . SO 

Ji 


where  the  c.  . are  rational  expressions  of  the  derivatives  of  p at  s. 

1 j • • • • 1 
J0 

Definition  4.1.  \!r  is  said  to  be  a normal  iteration  if  the  denominator  of 

each  c.  . is  a power  of  p' (s) . 

J0”*',Jm 

From  (4.1)  and  (4.2)  wc  have  that  both  Newton  iteration  and  secant  itera 
tion  are  normal.  In  fact,  most  commonly  used  iterations  v normal.  Wc 
prove  that  the  classical  one-point  inverse  interpolatory  iterations  \!*  (see 
Traub  [1964,  Section  5.1];  in  particular,  ^ *-s  Newton  iteration)  are 
normal.  Let  q denote  the  inverse  function  to  p and  v^  ^ = p(w^  ) . Then 


s = q (0)  = q(v(i))  - q'  (v(i))v(i)  + q"  (v(i) ) (v(  0 ) 2 + ... 


n"  m.  »v.1.'-'-  * »'-' 1 .-!  '.f'.'U“i^ .......  u, 


»o«- 
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By  definition  of  * , 
P 


= *(w(i))  + 'L  q(j)(v(i))(p(w(i)))j, 

P J • 


and 


j*P 


(i+1) . ^ (w<i>> <,<„<»>> j. 

L-‘  J* 
jsP 


Note  that 


p(w(i))  = P'(s)e(i)  + | p"(s)  (e(i))2  + ... 


and  that  q^(v^)  is  a rational  expression  of  p^(w^)  for  k=l,...,j  and 
has  the  denominator  (p1 (w(l)))j.  Expanding  the  p(k)(w(i))  around  s shows 


that  i|r  is  a normal  iteration. 


Definition  4.2.  For  a normal  iteration  i|r  defined  by  (4.3)  and  satisfying  (4.4), 
we  define  the  order  o of  by 


p = sup{r|rnvfl  £ jQrm  + 1 + . . . + Jm  for  all  (jQ,...,jm) 

such  that  c 


m-1 


jn>  • • • > J. 


. f 0 for  some  polynomial  p.} 
m 


By  (4.1),  it  is  easy  to  check  that  the  Newton  iteration  has  order  2. 

In  general,  it  can  be  shown  that  the  one- point  inverse  interpolator  iteration 

has  order  p.  Consider  now  the  secant  iteration.  By  (4.2),  the  order  of 
P 

the  iteration  is  given  by 

o - sup{r|r2  £ jr  + l for  all  j ,&  s l}, 

Q 

which  is  equivalent  to  o = sup{r|r  s r + l].  Hence  p is  the  positive  root 
of  r2  = r+1,  i.e,,  o = (j)  = (l+\^5)/2. 


**& 
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o.  REGULAR  PROBLEMS:  NORMAL  ITERATIONS  ON  POWER  SERIES 

We  show  how  normal  numerical  iterations  with  order  greater  than  unity  can 
always  compute  an  expansion  of  an  algebraic  function  for  a regular  problem. 
The  main  result  is  Theorem  5.1.  As  a corollary  of  this  theorem  we  show  that 
a Newton-like  iteration  always  "converges  quadratical ly" . Wo  also  show  the 
convergence  of  a secant-like  iteration.  We  end  the  section  with  an  example 
of  a convergent  first  order  iteration. 

We  begin  with  some  definitions.  Recall  that  a meromorphic  series  is  a 
power  series  with  a finite  number  of  negative  powers. 

Definition  5.1.  Given  a meromorphic  series  A(z)  and  a real  number  a,  then 
by  the  notation 

13(z)  = A(z)  (mod  z°) 


we  mean  B(z)  is  a finite  series  cons  is  tin,*;  of  all  terms  of  A(z)  of  degree  < a. 

Let  ft  be  a normal  numerical  iteration.  Let  the  numbers  w^,...,w^  m 
in  (A. .3),  the  defining  relation  for  ft,  be  replaced  by  meromorphic  '-cries 
(z) , . . . ,W^  m)  & . Then  the  iterate  W^^'"  (z)  defined  by 


«(i-"V.>  - 


is  in  general  a meromo* pa ic  series,  provided  that  it  is  vol  1-def in<  u.  Lot 
L',^  (z)  = w^1'^  (•'!  - S (z't  denote  the  error  <»:  the  i ti:  iterate. 


Definition  5.2.  Wo  say  an  iteration  on  meromorphic 


ri  es 


cos 


COllVi  ; 


holds  for  atroaorphic  series.  The  constant  tew  of  P'(S(z),z)  is  P*(sn,<0,  which 
is  non-zero  by  condition  (i) . Thus  by  conditions  (ii)  and  (iii),  (5.2) 


implies  diet 


ord(E(mfl))  * min(j0p®  + J^o""1  + ..  . * jm) 


where  the  minimum 


is  taken  over  all  the  such  that  C.  .is 

« m Ja»  > J 

J0  Jm 


non- rero  for  some  PtH,;:).  By  the' definition  of  p in  Section  4 , we  have 


> ip’*'. 


By  induction*  (5.1)  can  he  established  for  i^nf 1 ,mf2 , . . . , using  similar  argu- 
ments. The  convergence  of  the  iteration  follows  immediately  from  Definition 


Remark  5.1.  Thus  well-defined  normal  iterations  on  regular  problems  always 
converge.  This  behavior  is  strikingly  different  from  the  behavior  of  these 
iterations  on  numerical  polynomials  where  only  local  convergence  is  assured 
unless  Strong  conditions  are  imposed..  Note  that  the  expaf  >n  S(z)  may  cOn* 
verge  in  only  a small  disk  around  the  origin!  we  shall  not  pursue  'the  domain 
qf  convergence  hero.  ■ 


Remark  5.2.  (5.1)  shows  that  W is  a power  series  with  non-negative  powers 

only  rather  than  a meromorphic  series.  Until  this  fact  was  established  it 
was  necessary  to  work  over  the  field  of  mesomorphic  series,  ■ 


Remark  5 . 3.  Observe  that  we  do  not  define  order  for  power  series’,  valued 
iteration  but  only  for  normal  numerical  iterations. 
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Remark  5.4.  Note  that  in  Theorem  5.1  we  need  not  assume  that  An(0)  f 0. 
This  fact  will  be  used  in  the  proof  of  Theorem  6.1. 


We  apply  Theorem  5.1  to  two  specific  iterations.  We  begin  with  a 
Newton-like  iteration,  which  is  defined  by  (5.3)  below.  This  iteration  is  ob- 
tained from  the  numerical  Newton  iteration.  In  the  power  series  setting  we 
hesitate  to  call  it  Newton  iteration,  since  Newton  [1670]  actually  used  a dif- 
ferent method  for  computing  the  expansion.  His  method  computes  one  coeffici- 
ent per  iteration  and  in  general  is  not  as  efficient  as  the  Newton* like 
iteration  defined  below.  We  will  discuss  the  Newton- like  iteration  in  some 

detail  since  we  anticipate  it  will  be  one  of  the  most  commonly  used  iterations 

_ 5P 

in  practice.  Here  and  elsewhere  we  use  the  notation  P'  (W,z)  = -^<]W,z).  Re- 
call that  the  numerical  Newton  iteration  is  a normal  iteration  of  order  2. 
From  Theorem  5.1  we  have 


it 

Corollary  5.1.  If 

(i)  P(sQ,0)  = 0 and  P' (sQ,0)  0, 

(ii)  W^  = sQ,  i.e.,  ord(E^)  5 1, 

then  the  iterates  W^  generated  by  the  Newton- like  iteration, 


(5.3) 


W(i+1)(z)  ,W(i)(z) 


(i) 


(z) 


P'(W(l)(z),z) 


2 i+1 

(mod  z ) 


are  well-defined  and  satisfy 
(5.4)  ord(E(l))  ^ 2 1 


for  i=0,l,2,...,  and  hence  the  iteration  converges. 

& 

A result  similar  to  Corollary  5.1  has  been  proven  independently  by 
Professor  J.  Lipson  (Private  Communication). 


l.’roof. . He  need  only  show  that:  the  iterations  Vt'x/ (?'  < are  all  vdl-defined. 
This  holds  since  for  all  i the  constant  term  in  P*  («  ,z)  is  P*  (s0>0), 

which  is  non- zero.  ® 


Remark  5.5.  If  we  define  the  valuation  of  a power  series  A(z)  to  be  b 


•ord  (A.) 


where  b is  any  positive  constant,  then  Corollary  5.1  follows  from  a known 
theorem  in  valuation  theory  (see  Bachman  [1964,  Ch.  II,  Theorem  4 .2]) . 


Tt  is  easy  to  show  that  if  S(z)  is  a polynomial  of  degree  q,  then  itera- 
tion (5.3)  will  compute  it  in  9j  + * iterations.  By  a slight  modifica- 

tion of  the  hypotheses  of  Corollary  5.1  we.  can  replace  the  inequality  (5.4) 
by  equality. 


Corollary  5.2.  If 


(i)  P(sQ,0)  » 0,  P'(8fl,0)  0,  P"(s0,0)  0, 

(ii)  W(0)  = sQ,  ord  (E<0))  = 1, 


then  the  iterates  generated  by  the  Newton-like  iteration  satisfy  ord(Ev  ) = 


Corollaries  5.1  and  5.2  can  easily  be  generalized  to  any  ono-point  in- 
verse interpolatory  iteration  . 

As  our  second  example  we  consider  a secant- like  iteration.  One  has  to 
be  somewhat  careful  in  defining  this  iteration.  A straightforward  approach 


would  generate  iterates  by 

(i+1)  _ „(t)  w(0-w(i~D 

(5.5)  W = H 7n  HTir 

P(W'  ;)-.>(tTl  ') 


P(W^)  (mod  z^  ) 


where  <})  = (l-k/5)/2.  Then  becomes  undefined  when 


w<« 


happens  when  there  is  a "large"  gap  between  the  degrees  of  two  consecutive 
terms,  in  the  expansion  which  we  want  to  compute.  A solution  to  the  problem 
is  given  in  the  following  Corollary.  The  idea  is  to  use  a perturbed  in 
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ord(E^ : ^ ord(E(l)E(0' ) = ord(E(1))  + ord(E^) 

s min(ord(E(1)),F3)  + ord(p/0))  > + F2  = F^. 

By  induction,  one  can  similarly  prove  that  for  1=2,3,...,  is  well- 

defined  and  ord(E(x))  ? Fi+0.  ■ 

Results  similar  to  Corollary  5.3  hold  for  other  suitably  modified  itera- 
tions with  memory,  (i.c.,  iterations  with  m > 0 in  (4.3)). 

So  far  we  have  only  dealt  with  iterations  of  order  > 1.  Wc  now  consider 
an  iteration  with  order  one.  i».  fine 


„(i+D  - „o>  . p(»(1)) 
p-(w(0)> 


for  i=0,l,2,...  . Then 

,,  7)  AM)  m At)  . p,(a)e(l)  +lpM(s)e(i)  + - 

p'  (s')  + p"(s)e(0)  + .. . 


Bli2l  oC0)  (i> 

p'(s)  c 


E^Sl  (e^1))2  + 

2p'(s)  ( 


e^.^W 


i+ 

j>0,  tal 


where  the  c.  are  rational  expressions  wuose  denominators  are  powers  of 
J > *' 

p' (s)  . This  implies  that  the  iteration  is  normal  and  has  order  p = 1..  We 
may  use  the  iteration  on  power  series  and  obtain  the  following  theorem  which 
is  an  easy  consequence  of  (5.7): 

Theorem  5.2.  If 


(i)  P(s0,0)  = 0,  P'(s0,0>  4 o. 

(ii)  W(0)  = sn. 


I 


11 


6 THE  GENERAL  PROBLEM:  NEKTON  POLYGON  PROCESS 


Recall  that  our  general  problem  is  to  compute  the  value  of  d and  the  co- 
efficients s0'sl****  of  an  expansion 


of  the  algebraic  function  corres.x.nding  to  a given 


P (W,z)  = A (z)W  + ...  + An(z). 


In  this  section,  we  show  tnat  the  general  problem  can  be  reduced  to  a regular 
probl.-m  by  transforming  P(W,z)  to  some  P(W,z).  The  regular  problem  can  then 
be  :olved  by  normal  iterations,  as  described  in  Section  5. 

Since  P(s0,0)  = 0,  can  be  obtained  by  finding  a zero  of  the  numerical 
polynomial  P(W,0).  In  this  section  we  assume  that  finding  a zero  of  a numer- 
ical. polynomial  is  a primitive  operation.  (This  assumption  will  be  removed 
in  the  next  section  bv  carrving  the  zeros  .symbolically.)  If  P'  (Sq,0)  f-  0, 

wo.  have  a regu'.n  problem  solvable  by  a normal  iteration.  Hence  ve  assume 
that  P' (Sq,0)  ” 0.  i hen  is  a i*»ltiplc  v to  of  the  numerical  polyn.-”iaL 
P(W,0)  and  there  Is  more  than  one  < .".pension  of  tin  algebrai*  function  start- 
ing with  Sq.  Wo.  wool  ■ not  expect  an  iteration  starting  with  W^  - sn  to 
converge  since  the  iteration  would  not  "know'"  to  which  expansion  it  should 
converge.  Intuitively  the  convert',  nee  of  an  iteration  reqi ires  tint  it 
start  with  an  initial  segment  a unique  expansjoi  . 'Ibis  so  <..*  t(  that  v» 
find  an  initial  segment  of  a unique  expau:  ion  f rting  with  . fin*  existence 
of  the  segment  is  guaranteed  only  if  no  fa  expansions  coincid..  , i.e.,  o 


I.e.,  zeros  of  a polynomia1  can  he  computed  to  any  prespecified  prersi  a. 
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discriminant  D(z)  of  P(W,z)  with  respect  to  W is  not  identically  equal  to  zero. 
Therefore,  in  this  section  we  shall  assume  that 

D(z)  i-  0. 

The  assumption  holds  when  P(VJ,z)  is  irreducible  or  simply  when  P(W,z)  is 
square- free  (Walker  [1950,  Theorem  3.5]).  Hence  we  can  make  this  condi- 
tion hold  by  using  factorization  or  square-free  decomposition  algorithms 
but  do  not  pursue  this  here. 

A classical  method  for  finding  an  initial  segment  of  a unique  expansion 
uses  a geometric  aide  known  as  the  Newton  Polygon,  which  provides  a conveni- 
ent tool  for  analyzing  a set  of  inequalities.  (Some  authors  refer  to  Puiseux's 
Theorem  because  of  the  work  of  Puiseux  [1850]  but  clearly  the  idea  originated 
with  Newton  [1670,  p.  50].)  The  method  has  not  been  subject  to  algorithmic 
analysis . 

We  state  the  Newtcn  Polygon  Process  adapting,  with  some  modifications, 
the  description  in  Walker  [1950].  In  Theorem  6.1  we  show  that  the  Newton 
Polygon  Process  transforms  the  general  problem  to  a regular  problem.  Theorem 
6.1  also  gives  the  connection  between  the  number  of  identical  terms  in  at 
least  two  expansions  and  the  number  of  Newton  Polygon  stages.  Theorem  6.2 
gives  an  a priori  bound  on  the  number  of  stages  which  differs  by  at  most  a 
factor  of  two  from  the  optimal  bound.  Example  6.1  shows  that  in  general 
F(W,z)  must  be  transformed  to  a new  polynomial  P(W,z) ; it  is  not  enough  to 
compute  an  initial  segment  of  a uniqu-*  expansion  and  use  it  as  the  initial 
iterate  for  a normal  iteration  on  the  original  polynomial  P(W,z). 

In  the  following  algorithm,  let  A.  (z)  be  the  coefficient  of  W*  in 

Q\  . ’ 
i,k 

i,kz  be  the  lowest  degree  term  in  A (z) 

' .k 


Pk(W,z) . If  A.  k(z)  f 0,  let  a 


Newton  Polygon  Process 


Nl.  k «-  1,  Pk(W,z)  - P(W,z). 


N2.  Plot  the  points  f.  , = (i,<v.  . ) on  the  xy  plane  for  i such  that 

l)K  1,K 

A.  . (z)  0.  .loin  fn  to  f . with  a convex  polygon  arc  each 

l)K  u jK  n ,R  ' 

of  whose  vertices  is  an  f.  , and  such  that  no  f.  , lies  below  any 

i,k  L,lc 

line  extending  an  arc  segment. 

N3.  If  k a 1,  choose  any  «>«. tnt  y + v^x  * of  the  arc.  If  k **  1, 

choose  a segment  with  Vk>  0.  (Such  a segment  always  ■ s t s . ) 

Let  g,  denote  the  set  of  indices  i for  which  f.  . lies  on  the 
k i , k 

chosen  segment.  Solve  the  polynomial  equation 


a . , x =0. 
i,k 


Let  be  any  of  the  non- zero  roots.  (Such  a non- zero  solutic 
always  exists.) 

N4.  If  e is  a simple  zero,  go  to  Nb;  else  go  to  N5. 

- ft  y 

N5.  \+1(W,z)  *”  Z *^z  (Wfck),z),  k - k+1 . Go  t:o  N2 . 

N6.  t k.  (Hence  t represents  the  number  of  stages  taken  by  the 
Newton  Polygon  Process.) 

-?  y 

P(W,z)  - z ' • P (z  ^,2), 

P(W,z)  - P(W,zd), 


where  d is  the  smallest  common  denominator  of  y \ . (•■  may 

1 1 1 

be  zero.  If  y^  = 0 we  assume  that  - j has  one  us  its  denominator.) 
Terminate  the  process.  I 


Lemma  6.1.  After  the  Newton  Polygon  Process  terminates,  the  following 
properties  hold: 


(i)  The  coefficients  of  P(W,z)  are  polynomials  in  z. 

(ii)  is  a simple  zero  of  the  numerical  polynomial  P(W,0) 


Proof . It  is  easy  to  verify  (i).  To  prove  (ii)  we  show  that 


P(M,0)  = ^ a. >tw  . 


L 


For  notational  convenience,  let  a.  „ = <7. , a.  _ s a.,  B ~ fl,  V = v,  g 

i,t  1 i,t  i t Tt  " bt 

and  let  g denote  the  set  complementary  to  g with  respect  to  {0,l,...,n3. 


= 8 


Let 


O' 


Pt(W,z)  = (anz  + Qn(z))w“  + ...  +(aoz  + Qq(z)), 


where  ord (Q . ) > > . . Then 
u 1 


n 


i V—  3 

P (W , z)  - a.w1  + a.z  -1  W ! + zLV~  . (z)W1. 

i 4-  J - 1 


I'rg 


i-0 


Since  0 = + 1 


/ 0 


; 1 l-‘*  -L  8,  Vj  €•  g, 


P(W,U)  = P(W,0)  = aWl. 

i^g 

Theorem  t> . 1 . After  the  Newton  Polygon  Process  terminates,  the  following 
properties  hold: 

(i)  The  general  problem  of  computing  an  expansion  S(z)  of  the  algebraic 
function  corresponding  to  P(W,z)  has  been  reduced  to  the  to  I lowing 
rcgi lar  problem:  Compute  the  expansion  §(2)  starting  from  c 


for  the  algebraic  function  corresponding  to  Pft% ,z^ . Then 


let 


S(z) 


t-1 

U 

i=l 


■t  Y.+ . . .+ 

(ii)  S(z)  is  the  unique  expansion  with  starting  segment  ) c.z 

<Ji=l  1 

-j  YJ+...+V. 

(iii)  There  is  more  than  one  expansion  which  starts  with  j c z 

1=1 

for  every  j < t.  That  is,  there  are  at  least  two  expansions 
which  coincide  in  their  first  t-1. terms. 


Proof.  By  Lemma  6,1,  we  conclude  that  the  problem  of  computing  S(z)  is  regular. 
(Note  that  the  leading  coefficient  of  P(W,z)  may  vanish  at  z = 0.  See  Remark 
5.4.)  (i)  follows  from  P(W,z)  = P(W,z^)  and 


P(W,z)  = z 


t>‘  V"-+Yi  . Yi+---+Yt„  I 

c.z  + Z W,Z  I . 

^ 1 * 


(ii)  and  (ii-  » hold  since  the  Newton  Polygon  Process  does  not  terminate 
until  cfc  is  a simple  zero. 

*■  Vj+...-bV^ 

Since  there  is  only  one  expansion  which  starts  with  ^ c^z  , 

i=l 

we  might  expect  that  if  th'.s  segment  is  taken  as  the  initial  iterate 

for  a normal  iteration  then  the  iteration  on  the  original  polynomial  P(W,z) 


rather  than  on  the  transformed  polynomial  P(W,z)  will  converge.  The  follow- 
ing example  shows  this  not  to  be  the  case;  in  general  we  must  use  the  trans- 
formed problem. 


Example  6.1.  This  problem  appears  in  Jung  [1923,  p.  29]  although  it  is  not 
used  to  illustrate  the  point  we  wish  to  make  here.  Let 

2 3 12  4 

P(W,z)  = W - (2+z+z  )W  + 1 + z + -z  + z . 

The  two  expansions  are 

1 3/2  1 3/0 

S^(z)  =*  1 + ~z  + z + ...,  S2(z)  = 1 + ^z  -7.  + ...  . 

Suppose  that  we  want  to  compute  S^(z)  by  the  Newton-li!.e  iteration.  If  we 
take  W(0)  * 1 + jz  + z^2  in 


(i) 


„<i«> . „<»  . v>  , 


P'  (w(i)  ,z) 


we  fina  s 1 + tz  - 7Z^2  + ...  . differs  from  S,  even  in  the  co- 

2 4 1 

efficient  of  z^\  Though  there  is  only  one  expansion  starting  with  W^»  namely, 
Sj,  the  Newton- like  iteration  starting  from  does  not  converge  to  S^.  ■ 


We  illustrate  the  Newton  Polygon  transformation,  the  transformations  of 
Section  2 and  the  iterative  process  with  another  problem  in  Jung  [1923,  p.  31]. 


Example  6.2.  Find  all  the  expansions  of  the  algebraic  function  corresponding 
3 2 

to  P(W,z)  ■ -W  + zW  + z around  zQ  3 The  first  transformation  of  Section 

2 3 

2 converts  P(W,z)  to  -z  W + zW  + 1,  which  is  then  converted  by  another 

3 

transformation  to  -W  + zW  + z.  The  Newton  Polygon  Process  yields 


ypiuLuw  mums? 
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t = J,  0j  * 1,  Yj  **  1/3*  Cj  * i,  d « 3 and  P(W,z)  * -If*  + *H+  1.  Take 
W<°>  =»  1.  Then  the  Newton-like  iteration  (5.3)  applied  to  P(W,z)  gives 


Thus 


W(1)  = 1 + z/3,  W(2)  » 1 + z/3  - z3/81. 

S(z)  » z1/3S(z1/3)  = zl/3  + z2^3/ 3 - z4/3/81  + ...  . 


Let  T(z)  * S(z)/z  * z 2'^3  -f  z ^3/3  - z*/3/81  + ...  . Then  an  expansion  of 
the  given  problem  is 


T(- 


i)  . ,2/3  + 11/3  .1-1/3  + _ . 


81 


The  other  two  expansions  are 


0z 


2/3  £ 1/3 . 8.  -1/3 

* 7 *•  81  * • • • • f 


.2  2/3  ^ e 1/3  0z  -1/3  , 

0 z +-JZ  "81Z  +...» 


where  0 is  the  primitive  third  root  of  unity. 


The  following  theorem  gives  an  a priori  bound  on  the  number  t of  stages 
in  the  Newton  Polygon  Process  which  differs  by  at  most  a factor  of  two  from 
the  optimal  bound. 

Theorem  6.2. 

(6.2)  t < ord(D)  + 1 

Furthermore  for  all  t there  exist  problems  for  which  t ■ ord(D). 

Proof,  The  theorem  is  trivial  if  t = 1.  We  assume  that  t ^ 2.  Then  by 
(iii)  of  Theorem  6.1,  there  are  at  least  two  series  expansions  S,  ami 
which  agree  in  the  first  t - 1 «cn-zero  terms.  Write 


•Sk*'--^* * »^^a-^r--;o<»W!|»t.gt.L^.«V  uW't?g.jK{« -.ft-.i^/...'  \M.iy.  »»w.»  ; •■  ^»vV 
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S,  = s -z  , 
1 — l,a.  * 

i=l  1 


s2=  -.s2,b;z  > 

i=l  1 


where  the  {a^},  {b^}  are  strictly  increasing  non-negative  integer  sequences 

such  that  none  of  the  s.  , s vanish  and  s = s , , a /d,  = b./d_ 

2 ,£k  2,tu  l 1 i 2 

for  i=l,...,t-l.  Without  loss  of  generality,  assume  d^  £ d^.  Note  that  the 
cycle  which  contains  has  the  series: 


■H  d, 

■i.j  * ,.*i,«fh  '*  > J"0-' 

i=l 


. »d,~l  , 


and  the  cycle  which  contains  has  the  series: 


S2,J 


i=l 


b. 

jbi  T2 

s 2 ^ ’z  ’ » • • • > ^2~ ^ * 


where 


2Ui£l 

dl  d2 

?,  ■ e and  i,  “ e 


Note  that  we  do  not  rule  out  the  possibility  that  S.,  and  S£  are  in  the  same 
cycle  and  that  therefore  the  cycles  (S^}  and  are  identical.  Si 


nice 


2tk£T 


ja. 


d.  jai  jb 


p 1 , 1 

h “ e 


^2  1 = 6 2 


2tt,CT., 

T~jbi 


and  a./d,  - for  W t-1.  S^j  a"*!  *2>J  to  the  first  t-1  terms 


for  j=0, . . . ,d1-l. 


Hence , 


U4vy  y c^gawwg^gwgg  ■ wjffg.vi'.'V'w.^ 


'wv  »■ , .«.  ^wjih^li  ^gy  ji 


I- 
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t 

? 

I 

< 

!? 

ijf 

< 

I 


OTd(si,j-f'2,j! 


• bt-i+1  A'.I 

" d2  ' dl ' ’ d2 


]>t 


Then 


dr1 

»(*>  * n (s,  ,(z)-s„  ,(z)). 

j=fl  i»J  z»-> 


ord(D)  £ ord(V) 


t- 1 1 

sdi(4-i  + -H 

1 dj  d2 

2 Vi  + '• 

Since  the  {a^}  is  a strictly  increasing  non-negative  integer  sequence, 
at^  £ t-2.  Thus,  ord(D)  > t- 1 which  establishes  (6.2).  Let 

t 

Sj(x)  •-  ^ z'\  S2(a)  = S^z)  - zt 

j=0 

and 


P(W,z)  = (W-S^zJXw-S-Cz)). 

A / 

2 

By  Theorem  6.1,  the  Newton  Polygon  Process  has  t stages.  ord(D)  = ord^Sj-Sg)  ) 
which  completes  the  proof.  ■ 

Theorem  6.2  gives  a c*'m:  utable  a priori  bound  but  requires  the  computa- 
tion of  ord(D).  A very  cheap  bound  is  given  by 


= 2t 


Corollary  6.1. 


t < m(2n- 1)  + 1 


where  m **  max  (deg  A.). 
nsi<n 


Proof,  d(z)  is  a determinant  of  order  2n-l  whose  elements  are  polynomials 
of  maximal  degree  m.  Hmcc  D(z)  is  a polynomial  of  degree  at  most  m(2n-l). 
Since  D(z)  cannot:  vanish  identically,  ord(D)  m(2n-l).  I 


-SKf 
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7.  A SYMBOLIC  MODE  OF  COMPUTATION 

The  Newton  Polygon  Process  involves  computing  roots  of  polynomial  equa- 
tions (6.1).  Instead  cf  actually  solving  the  equations,  in  this  section  we 
carry  the  roots  symbolically  through  their  minimum  polynomials.  We  assume 
that  the  underlying  field  A is  one  where  exact  arithmetic  can  be  performed  such 
as  a finite  field  or  the  field  Q of  rational  numbers.  Then  the  expansions  can 
be  computed  symbolically  with  exact  arithmetic.  The  following  example,  where 
A is  taken  to  be  Q,  will  illustrate  the  idea. 

Example  7.1. 

P(W.z)  = W3  + (z+z2)W2  - 2z2W  - 2z3. 

We  snail  compute  an  expansion  of  the  algebraic  function  corresponding  to 

P(W,z),  using  exact  rational  arithmetic.  The  first  stage  of  the  Newton 

3 2 

Polygon  Process  yields  y^  = 1,  3^  = 3 and  c^+c^-  2c ^ - 2=0.  Since 
c3  + c2  - 2c1  - 2 = (c2-2) (c^+1) , cx  = • £,  -J5  or  -1.  Suppose  that  we  are 
interested  in  the  expansion  starting  with  Ji  or  -Ji,  Instead  of  using  an 

approximation  to  JZ  or  we  carry  Cj  symbolically  through  its  minimal 

2 

polynomial  M^(x)  = x - 2.  That  is, 

(7.1)  c2  - 2 » 0. 

Since  the  equation  has  only  simple  zeros,  the  Newton  Polygon  Process  termi- 
lates  with  t = 1,  and 

P(W,z)  = z"3P(zW,z) 

= W3  + (l+z)W2  - 2W  - 2. 

We  use  the  Newton-like  iteration  (5.3)  to  compute  S(z)  such  that  P(S(z),z')  0. 
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Let  W^(z)  * Cj-  Then 


c*+(l+z)c  2-2c  -2 

WU;(z)  = c-  - -i-r i 

Zcf2(\+z)cx-2 


2 

(mod  z ) 


Using  (7.1),  we  obtain 


W^(z)  = ~ -|z. 


Similarly  all  coefficients  of  z-^  in  W^(z)  can  be  represented  as  linear 
polynomials  in  with  rational  coefficients.  By  (il)  of  Theorem  6.1,  a 


solution  to  the  given  problem  is 


1 2 


S(z)  “ TS(2>  3 CjZ  - -z  + ...» 

which  represents  both  the  numerical  expansions  starting  with  Jlz  and  -Jlz.  ® 

In  gene'al,  when  the  Newton  Polygon  Process  is  performed,  c^>  k=l,...,t, 
can  be  carried  symbolically  through  its  minimum  polynomial  M^(x)  over 
Q(Cj,..,,c^j).  Then  all  the  coefficients  of  the  expansion  S(z)  are  in  the* 
extension  field  Q(c1> . . . ,ct> . To  simplify  the  computation,  one  can  compute 
from  M^(x)  the  minimum  polynomial  M(x)  for  c,  where  c is  a primitive  element 
of  the  extension  field  Q(c^, . . . ,0^) , i.e.,  Q(c)  = Q(c . »c^) . Then  the 
coefficient^  of  the  expansion  S(z)  can  all  be  represented  by  polynomials  of 
the  form  wi'ere  h - dog  M and  € Q.  S(z)  can  be  computed  entirely 

with  exact^arrthmetic  . Furthermore,  S(z)  give  a simultaneous  representation 
•of  h numerical  expansions;  S(z)  can  be  used  to  produce  h numerical 
expansions  by  substituting  zeros  of  M(x)  for  c in  the  coefficients  of  S(z'). 
(This  implies  that  h s n.) 


JP.PUJ.',.  .J'U.l,..  1 . '■ 1 -V  ..  —u,  l.-J. 
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8.  ASYMPTOTIC  COST  ANALYSIS 


In  this  section  we  analyze  the  cost  of  computing  the  first  N terms 
(including  zero  terms)  of  an  expansion  for  large  N.  Since  the  Newton 
Polygon  Process  is  independent  of  N.  by  Theorem  6.1  we  can  without  loss  of 
generality  assume  the  problem  is  regular.  Furthermore,  since  the  asymptotic  re- 
sults will  be  the  same  for  any  normal  iteration  with  order  greater  than  one, 
we  shall  assume  that  the  iteration  (5.3)  is  used.  Our  cost  measure  is  the 
number  of  operations  used  over  the  field  A.  If  we  carry  zeros  symbolically  as 
described  in  Section  7,  then  we  work  over  an  extension  field  A(c)  rather 

than  A.  If  the  minimum  polynomial  for  c is  of  degree  h,  then  operations 

2 

in  A(c)  are  more  expensive  than  in  A by  a factor  of  0(h)  or  0(h  ). 

Since  h is  independent  of  N,  in  our  analysis  we  shall  not  be  concerned  with 
whether  or  not  zeros  of  polynomials  are  carried  symbolically. 

Let  M(j)  be  the  number  of  operations  needed  to  multiply  two  jth  degree 
polynomials  over  the  field  A.  Assume  that  M(j)  satisfies  the  following  mild 
condition:  there  are  a,  (3  € (0,1)  such  that 


(8.1)  M(  r<yj! ) ^ pM(Pfl) 

for  all  sufficiently  large  j.  Observe  that  W^(z)  is  a polynomial  of  degree 
at  most  2*  - 1,  and  that  the  computing  \i^+L\z)  by  (5.3)  takes  0(nM(2i-l)) 
operations.  Hence  the  total  cost  of  computing  N terms  in  the  expansion  is 
0(n(M(N)  + M(fN/2l)  + M(rN/4l)  + ...)),  which  is  0(nM(N))  by  condition  (8.1).  (See 
Brent  and  Kung  [ 1976,  Lemma  1.1].)  We  summarize  the  result  of  this  section 


in  the  following 


i 


yi/y^y-r'3!ir.'a  w 


9.  EXAMPLES 


We  choose  as  our  examples  calculation  of  the  Legendre  polynomials 
through  their  generating  function,  solution  of  an  equation  with  transcen- 
dental coefficients,  and  calculation  of  the  expansion  of  a complete  ellip- 
tic integral.  Although  the  first  two  examples  are  not  covered  by  the  theory 
of  this  paper,  they  are  covered  by  easy  extensions  of  our  results.  Examples 
9.1  and  9.3  are  illustrations  of  the  many  applications  of  algebraic  function 
expansions. 

We  use  the  Newton-like  iteration  (5.3)  in  all  three  examples  with  the 
notation: 


P.  a P(W(l)  (z)  ,z, , P^h  P*  (W(l)  (z)  ,z)  = -~y  (WU)  ,zS  , ^ = 

Within  each  iteration  step  we  exhibit  enough  terms  so  that 
computed  to  2*-l  terms. 


(z)  can  be 


Example  9.1.  Legendre  Polynomials 

The  generating  function  for  Legendre  polynomials, 

1 “ 

(l-2tz+z2)  2 = L.Ctlz1 
i=0 

satisfies 

P(W,z,t)  = (l-2tz+z2)W2  - 1. 

Take  W(0)  =>  1.  Then 

PQ  = -2tz , P|j  * 2,  6Q  = -tz,  W(1)  =■  1+tz. 

P:  = (l-3t2)z2  + (2t-2t3)z3,  P^  = 2(l-tz) , 6j  -|(l-3t2)z2  + |(3t-5t3)z\ 

(2)  12  2 13  3 

W - 1 + tz  + j(3f-l)z  + 2 (5t  -3t)z  . 


Hence  the  first  four  Legendre  polynomials  are 


L (t)  = 1,  L (t)  * t,  L (t)  »-|(3t2-l)  and  L-(t)  =*  j(5t3-3t) 


B.  Neta,  a student  at  CMU,  computed  the  first  32  Legendre  polynomials 
by  this  iteration  using  MACSYMA. 


Example  9.2 


P(W,z)  = W + (z+l)W  + sin  z. 


3 5 7 

7,  Z 7,  ( 0) 

Note  that  sin  z = z-  -^-  + tt-'57+...  . Take  W 

- ' • J « / • 


0.  Then 


0 


7-,  Pq  = 1,  60  = z,  W 


i 


0) 


-Zi 


pi”>  w<2)  --z  + r- 


Example  9.3.  A Complete  Elliptic  Integral 
Define  the  integral,  by 


Let 


"2  2 2 ~ 1 /9 

£(t)  - j (1-t4  sin4  0)  1 4d0. 

*n 


P(W,z)  = (l-z)W2  - 1,  z = t2  sin2  0. 


Take  = I.  Then 


0 0 “ 5 0 2 ’ 1 2 * 


„ . 12  1 3 n,  _ 0 , 3 2 53 

_ ~ uz  * ri  2 " z»  'i  "sz  ” \(7‘  ’ 


w 


(2^i+f  + |Z; 


j — — — •> 


in"  * 


W^4  is  an  initial  segmo  t of  the  algebraic  function  S(z)  correspond  inn  to 


P(W,z).  Since 
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n 

f(t)  = J2S(t2  sin2  0)dfi, 

0 

f(t)  = T)0  + |l)1 1 2 + |ll2t4  + ^t6  + .... 

where 

TT 

T1  [*2  . 2i 

'I.  - J sm  9d9. 

1 0 

For  this  simple  example  the  result  can  be  obtained  directly  by 
expansion  but  this  cannot  of  course  be  done  in  general. 


a binomia 


i 


10.  EXTENSIONS 

Our  ala  la  tills  paper  has  been  to  show  that  algebraic  functions  fora 
Interesting  and  useful  domain  in  which  to  do  algorithmic  complexity 
analysis  and  to  exhibit  fast  algorithms  for  computing  any  expansion  of  an 
algebraic  function.  In  this  initial  paper  we  have  restricted  ourselves  to 
the  "pure"  case  of  algebraic  functions  where  P(W,z)  is  a polynomial  in  W 
with  polynomial  coefficients.  We  list  some  additional  problems  which  we 
hope  to  discuss  in  the  future.  For  a number  of  these  out  results  (especi- 

ally  on  regular  problems)  apply  with  minor  modifications;  others  will  require 
major  new  results.. 

1.  Let  W be  a scalar  variable  but  take  z to  be  a vector  variable.  Re- 
sults similar  to  those  in  Section  5 should  hold.  We  have  seen  this 
case  in  Example  9.1; 

2.  Let  the  coefficients  of  P,  A^(z) , be  power  series  (rather  than  poly- 
nomials) . Results  similar  to  those  in  Section  5 should  hold.  See 
Example  9.2, 

3.  Let  both  W and  z bo  vector  variables.  This  is  the  fully  multivari- 
ate case,  which;  except  for  regular  problems,  is  iri  general  very 
difficult, 

4 , The  domain  over  which  we  have  worked  is  hot  algebraically  , closed 

since  problems  with  polynomial  coefficients  lead  to  solutions  repre- 
sented by  fractiprial  power  series  . If  the  coefficientr,  are  frac- 
tional power  series,  the  domain  is  algebraically  closed  (Puiseux’s 
theorem,  see,  e.g. ; Lefschetz  [1953])  and  tills  is  therefore  a 
natural  setting.  The  Newton- like  iteration  is  still  valid  oh  frac- 
tional power  series  for  .regular.  problems. 


5.  The  field  A need  not  be  restricted  to  the  complex  number  field. 

It  is  of  particular  interest  to  extend  all  the  results  to  finite 
fields. 

6.  An  important  computational  model  is  the  "fully  symbolic"  ode  where 
the  coefficients  of  the  expansion  series  are  expressed  as  functions 
of  the  input  coefficients. 


7.  Perform  complexity  analysis  which  includes  the  cost  due  to  the 
"growth"  of  coefficients. 
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